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Abstract 

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for 
Dark Matter. For understanding the nature of WIMPs and identifying them among new 
particles produced at colliders (hopefully in the near future), determinations of their mass 
and couplings on nucleons from direct Dark Matter detection experiments are essential. 
Based on our model-independent method for determining the WIMP mass from experi- 
mental data, I present a way to also estimate the spin-independent (SI) WIMP-nucleon 
coupling by using measured recoil energies directly. This method is independent of the 
velocity distribution of halo WIMPs as well as (practically) of the as yet unknown WIMP 
mass. In a background-free environment, for a WIMP mass of ~ 100 GeV the SI WIMP- 
nucleon coupling could in principle be estimated with an uncertainty of ~ 15% by using 
2 (or 3) X 50 events from experiments. 



1 Introduction 



Astronomical observations and measurements indicate that more than 80% of all matter in the 
Universe is dark (i.e., interacts at most very weakly with electromagnetic radiation and ordinary 
matter) . The dominant component of this cosmological Dark Matter must be due to some yet to 
be discovered, non-baryonic particles. Weakly Interacting Massive Particles (WIMPs) x arising 
in several extensions of the Standard Model of electroweak interactions are one of the leading 
candidates for Dark Matter. WIMPs are stable particles with masses roughly between 10 GeV 
and a few TeV and interact with ordinary matter only weakly (for reviews, see Refs. [1, 2]). 

Currently, the most promising method to detect different WIMP candidates is the direct 
detection of the recoil energy deposited in a low-background underground detector by elastic 
scattering of ambient WIMPs off target nuclei [3, 4]. The recoil energy spectrum can be calcu- 
lated from an integral over the one-dimensional velocity distribution function of halo WIMPs, 
fi{v), where v is the absolute value of the WIMP velocity in the laboratory frame. In our earlier 
work [5] , we presented a way to reconstruct this one-dimensional velocity distribution function 
and to estimate its moments from the recoil spectrum as well as from measured recoil energies 
directly in direct Dark Matter detection experiments. Neither the WIMP-nucleus scattering 
cross section nor the local WIMP density is required in this analysis. 

However, the mass of halo WIMPs is needed for the reconstruction of the (moments of the) 
WIMP velocity distribution. Therefore, as the next step we developed a model-independent 
method based on the reconstruction of the moments of fi{v) for determining the WIMP mass 
by combining two sets of (future) experimental data with different target nuclei directly 
[6, 7]. To do so, one simply requires that the values of a given moment of fi{v) estimated 
by both experiments agree. This leads to a simple expression for determining m^, which can 
be solved analytically and each moment can be used. Moreover, by assuming that the ratio 
of the spin-independent (SI) scattering cross sections on protons and on neutrons is known, 
an additional expression for determining has be derived. By combining the estimators for 
different moments with each other and with the estimator derived by making the assumption 
about the ratio of the SI cross sections, one can yield the best estimate of the WIMP mass [7] . 
Here we found again that neither a prior knowledge about the WIMP-nucleus cross section nor 
that about the local WIMP density is required. 

Meanwhile, in the second method for the determination of the WIMP mass, the product 
of the local WIMP density times the SI WIMP-proton cross section, PoO'^pj appearing in the 
expression for the scattering spectrum cancels out when we use the identity of this product 
for two different targets. Hence, as will be shown in the paper, once the WIMP mass can be 
determined one could then use this information to estimate cr^p conversely. Remind that, in 
order to identify new particles produced at e.g., the Large Hardon Collider (LHC) to be indeed 
WIMPs detected by direct detection [8], estimates of or constraints on their mass and couplings 
on nucleons from direct detection experiments are essential. However, due to the degeneracy 
between po and u^p, for estimating the SI WIMP cross section by this method one has to make 
an assumption for the local WIMP density, which can so far be estimated with an uncertainty 
of a factor of ~ 2 [1, 2]. Nevertheless, our simulations show that, in spite of the large statistical 
uncertainty due to very few events, for a WIMP mass of ~ 100 GeV, cr^p could be estimated with 
an uncertainty of 30% by using 2 (or 3) x 50 events from experiments. This result is (much) 
better than our estimate of the local Dark Matter density. 

The remainder of this article is organized as follows. In Sec. 2 I discuss the possibility of 
constraining the WIMP mass and its coupling on nucleons from a single experiment. In Sec. 3 I 
present the method for estimating the spin-independent WIMP-nucleon coupling by combining 



two (or more) experiments. Some numerical results based on Monte Carlo simulations of future 

experiments will also be presented. In Sec. 4 the analysis will be extended to the case of spin- 
dependent (SD) WIMP-nucleon couplings. I conclude in Sec. 5. Some technical details for our 
analysis will be given in an appendix. 



2 Constraining the SI WIMP— nucleon coupling 

The basic expression for the differential event rate for elastic WIMP-nucleus scattering is given 
by [1]: 



dv . 



dQ 

Here R is the direct detection event rate, i.e., the number of events per unit time and unit mass 
of detector material, Q is the energy deposited in the detector, F{Q) is the elastic nuclear form 
factor, fi{v) is the one-dimensional velocity distribution function of the WlMPs impinging on 
the detector, v is the absolute value of the WIMP velocity in the laboratory frame. The constant 
coefficient A is defined as 

^-7^, (2) 

where po is the WIMP density near the Earth and (Tq is the total cross section ignoring the form 
factor suppression. The reduced mass mr,N is defined by 

mr,N = : , (3) 

where is the WIMP mass and ttin that of the target nucleus. Finally, fmin is the minimal 
incoming velocity of incident WIMPs that can deposit the energy Q in the detector: 

Vmin = a^/Q , (4) 
with the transformation constant 



and I'max is the maximal WIMP velocity in the Earth's reference frame, which is related to the 
escape velocity from our Galaxy at the position of the Solar system, Vesc ^ 600 km/s. 

The local WIMP density at the position of the Solar system, po, appearing in the expression 
(1) for the scattering event rate has conventionally been determined by means of the measure- 
ment of the rotation curve of our Galaxy. The currently most commonly used value for po is 
[1,2] 

Po ~ 0.3 GeV/cm^ . (6) 

However, as mentioned in the introduction, due to our location inside the Milky Way, it is more 
difficult to measure the accurate rotation curve of our own Galaxy than those of other galaxies; 
an uncertainty of a factor of ~ 2 has thus usually been adopted [1, 2]^: 

Po = 0.2 - 0.8 GeV/cm^ . (10) 



^Recently, some new techniques have been developed for determining po with a higher precision [9, 10, 11, 12, 
13]. These estimates give rather larger values for po; e.g.. Catena and Ullio gave [9] 

Po = 0.39 ± 0.03 GeV/cm^ , (7) 



On the other hand, in most theoretical models, the spin-independent WIMP interaction 
on miclcus with an atomic mass number A > 30 dominates over the spin-dependent (SD) 
interaction [1, 2]. Additionally, for the lightest supersymmetric neutralino, which is perhaps 
the best motivated WIMP candidate [1, 2, 15], and for all WlMPs which interact primarily 
through Higgs exchange, the SI scalar coupling is approximately the same on both protons p 
and neutrons n [16]. The "pointlike" cross section ao in Eq. (2) can thus be written as 

-4 



2 



and the SI WIMP cross section on protons (nucleons) can be given as 

4^=(^)<l/p^ (12) 

where /p(n) are the effective xxpp(nn) four-point couplings, A is the atomic mass number of the 
target nucleus, and mr,p is the reduced mass of the WIMP mass and the proton mass rup. 
Here the tiny mass difference between a proton and a neutron has been neglected. 

As mentioned in the introduction, in our earlier work it has been found that one could in 
principle determine from direct detection experiments with neither a prior knowledge of ao 
nor that of po [6, 7]. Conversely, I will show in this article that one could also estimate or at least 
constrain the WlMP-nucleon cross section from experimental data directly without knowing rriy^, 
but for this estimation an assumption about po is needed. 



2.1 Expression for estimating the SI WIMP— nucleon coupling 

Our analysis starts from the expression (1) for the event rate for the elastic WIMP-nucleus 
scattering directly. By using a time-averaged recoil spectrum, and assuming that no direc- 
tional information exists, the normalized one-dimensional velocity distribution function of halo 
WIMPs, has been solved from Eq. (1) analytically [5] and, consequently, its generalized 

moments can be estimated by [5, 7]^ 

rv(Qina.x) 

K)(t'(gmin),^^(gmax))= / fi{v) dv 

and Salucci et al. even gave [11] 

po = 0.43 ± 0.11 ± 0.10 GeV/cm^. (8) 

Moreover, instead of a spherical symmetric density profile assumed in Refs. [9, 11], in Refs. [10, 12, 13] the 
authors considered an axisymmetric density profile for a flattened Galactic Dark Matter halo [14] caused by the 
disk structure of the luminous baryonic component. It was found that the local density of such a non-spherical 
Dark Matter halo could be enhanced by 20% or larger [10, 12] and Pato et al. gave therefore [12] 

Po = 0.466 ± 0.033(stat) ± 0.077(syst) GeV/cm^ . (9) 

^Here we have implicitly assumed that Qmax is so large that terms involving —2Qmax ^^^r{Qnia.x)/F^{Qma,x) 
are negligible. Due to sizable contributions from large recoil energies [5], this is not necessarily true, especially 
for some not very high Qmax in the experimental reality, and/or heavy detector targets, and/or heavy WIMPs. 
Nevertheless, considering the large statistical uncertainties due to (very) few events in the highest energy ranges, 
this should practically be a good approximation. 
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(13) 



Here ^(Q) = ot\/Q, Q(min,max) the experimental minimal and maximal cut-off energies of the 
data set, respectively, 



\"^/ expt, Q=Qn 



(14) 



is an estimated value of the measured recoil spectrum {dR/dQ)cxpt {before normalized by an 
experimental exposure, S) at Q — Qmin, and /„((5min, Qmax) can be estimated through the sum: 
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where the sum runs over all events in the data set that satisfy Qa G [Qmm, Qmax] and A'tot is the 
number of such events. Note that, firstly, by using the second line of Eq. (13) (i'")(t'((5min), ^(<3max)) 
can be determined independently of the local WIMP density po, of the velocity distribution 
function of incident WIMPs, fi{v), as well as of the WIMP-nucleus cross section ctq. Secondly, 
^^(Qmin) and In{Qmm,Qma.x) are two key quantities for our analysis, which can be estimated ei- 
ther from a functional form of the recoil spectrum or from experimental data (i.e., the measured 
recoil energies) directly^. 

By substituting the second expression in Eq. (11) into Eq. (1), and using the fact that the 
integral over the onc-dimensional WIMP velocity distribution on the right-hand side of Eq. (1) 
is the minus-first moment of this distribution, which can be estimated by Eq. (13) with n — —1, 
we have 
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(16) 



Using the definition (5) of a, the squared SI WIMP coupling on protons (nucleons) can be 
expressed as [17] 



If ( I ^ 



mil 

F (Qmin) 



(17) 



Note that, firstly, the experimental exposure E appearing in the denominator relates the actual 
counting rate {dR/dQ)expt to the normalized rate in Eq. (1). Secondly, due to the neglect of the 
terms -2Q^/^r(Qmax)/-P'^(Qmax) and -2r(Qmax)/-P'^(Qmax) in the denominator and numerator 
of the expression (13) for {v"'), respectively, |/pp determined by Eq. (17) would be overestimated, 
since the contributions from the two neglected terms are negative and the former is much larger 
then the later. However, because |/pp estimated by Eq. (17) is inversely proportional to the local 
WIMP density, whose commonly used value would possibly be underestimated (see Eqs. (6) to 
(10)), one should therefore at least be able to give an upper bound on |/pp. Then, by using the 



•^AU formulae needed for estimating r((5min), In{Qmin,Qn 
appendix. 



{), and their statistical errors are given in the 



standard Gaussian error propagation, the statistical uncertainty on |/pp estimated by Eq. (17) 
can be given as 




2A/'mC0v(m^, I /Mm) 
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2.2 From a single experiment 

The expression (17) for estimating the (squared) SI WIMP-proton couphng depends on three 
quantities: r((5min)5 -^o^ and the WIMP mass m^. As argued in Ref. [7], from a single recoil 
spectrum one cannot estimate without making some assumptions about the velocity distri- 
bution fi{v). Hence, as a model-independent analysis, one could only express/constrain |/pp 
as a (hnear) function/interval of the WIMP mass on the couphng-mass plane by using Eq. (17) 
with a single experiment. Meanwhile, from Eqs. (1) and (2), it can be found that, due to the 
degeneracy between the local WIMP density po a-^d the WIMP-nuclcus cross section (Tq) 
cannot estimate both of them independently^. Thus, for using Eq. (17), the simplest way is 
making an assumption for the local WIMP density po- 

In Fig. 1 I show the simulated results for a ^^Ge target with 5,000 experiments based on 
the Monte Carlo method^. The theoretical predicted recoil spectrum for the shifted Maxwellian 
velocity distribution [1, 2, 5] with a Sun's orbital velocity in the Galactic frame Vq = 220 km/s, 
an Earth's velocity in the Galactic frame = 1.05 fo,^ and a maximal cut-off velocity of the 
velocity distribution function v^ax — 700 km/s, as well as the commonly used elastic nuclear 
form factor for the SI cross section [20, 1, 2]: 



have been used. The SI WIMP-proton cross section has been set as 10^^ pb. The commonly 
used value of po = 0.3 GeV/cm^ has been used for both predicting the recoil spectrum and 
analyzing generated events. The experimental maximal cut-off energy Qmax has been set as 
50 keV and the threshold energy has been assumed to be negligible. Each experiment contains 
an expected number of 50 total events; the actual event number is Poisson-distributcd around 
this expectation value. The mass of incident WIMPs has been chosen as 25 (dotted magenta), 
100 (dashed blue), and 300 (double-dashed black) GeV, respectively. 

As we can see here, the prefactor, i.e., the slope of the linear function \fp\'^{'m^), in Eq. (17) 
is obviously underestimated. For the case of an input WIMP mass of 300 GeV, the theoretical 
value of |/pP (the filled green square) is even outside the la statistical uncertainty interval. 
This is because that the experimental maximal cut-off energy has been set as only 50 keV here. 
Remind that it is usually assumed that the WIMP flux on the Earth is negligible at velocities 

^In contrast, as I will show in Sec. 4, the ratios between different WIMP-nucleon couplings/cross sections can 
be determined without knowing the mass and the local WIMP density [18, 17, 19]. 

^Note that, rather than the mean values, in this article we give always the median values of the reconstructed 
results from the simulated experiments. 

^The time dependence of the Earth's velocity in the Galactic frame [1, 2] has been ignored. 
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Figure 1: The squared SI WIMP-proton couplings |/pp estimated by Eq. (17) and the lower 
and upper bounds of their la statistical uncertainties as functions of the WIMP mass for a ^^Ge 
target. The theoretical predicted recoil spectrum for the shifted Maxwellian velocity distribution 
with vq = 220 km/s, Ve = 1.05 vq, and fmax = 700 km/s as well as the commonly used elastic 
nuclear form factor for the SI cross section given in Eq. (20) have been used. The SI WIMP- 
proton cross section has been set as 10~® pb. The experimental maximal cut-off energy Qmax has 
been set as 50 keV and the threshold energy has been assumed to be negligible. Each experiment 
contains 50 total events on average. The mass of incident WIMPs has been chosen as 25 (dotted 
magenta), 100 (dashed blue), and 300 (double-dashed black) GeV, respectively. The filled green 
squares indicate the input WIMP masses and the theoretical values of |/pp. See the text for 
further details. 



exceeding the maximal velocity fmax- This leads thus to a kinematic maximum of the recoil 
energy 

max, kin n • \^^) 

For a WIMP mass of 100 (300) GeV, this kinematic maximum for a Ge target is 264 (504) keV. 
Hence, Jq in the prefactor of the linear function \ fp\'^{m^) given in Eq. (17) has been (strongly) 
underestimated. In Fig. 2 we increase therefore the maximal cut-off energy Qmax to 100 keV. 
It can be seen clearly that, by extending the detector sensitivity to higher energy ranges, the 
underestimated Iq and thereby the prefactor of the linear function |/pP(^x) '^^^ corrected 
significantly''. 



'''Remind that, since we neglected the term — 2(5max?'((5max)/^'^(Qmax) in the second bracket in Eq. (17), 
which contributes negatively, all results shown in this paper are somehow overestimated. 



^'^Ge, Q^ax < 1 00 keV, 50 events, a J^' = 1 pb 
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Figure 2: As in Fig. 1, except that the maximal cut-off energy Qmax has been increased to 
100 IceV. 



On the other hand, by substituting Eq. (17) into Eq. (12), one can express the SI WIMP- 
proton (nucleon) cross section as a function of the WIMP mass, a^p{m^), on the cross section 
versus WIMP mass plane. In Figs. 3 I show the simulated results for a Ge target with an 
input WIMP mass of 100 GeV. The experimental minimal cut-off energy Qmm has been set 
as (upper) and 5 (lower) keV. As a comparison I show also four extra curves drawn conven- 
tionally by using the shifted Maxwellian velocity distribution with four Sun's orbital velocities: 
vo = 180 km/s (dash-double-dotted orange), vq = 200 km/s (dash-dotted cyan), vq = 220 km/s 
(dotted magenta), vq = 240 km/s (double-dotted black), and the form factor given in Eq. (20). 

As shown here, two results analyzed by Eq. (17) and by the conventional method with an 
assumed halo model are compatible with each other in the mass and cross section ranges around 
and higher than the input values, whereas in the low WIMP mass range, these two curves show 
a significant incompatibility. Hence, by comparing results from these two analyses, one could 
in principle - for the first step with only one experiment observing positive signals - give the 
lower bounds of the WIMP mass and its cross section on protons (nucleons) (m^ > 40 GeV and 
(j^p ^ 7 X 10^^ pb from the upper frame of Figs. 3 in our simulation) from a single experiment. 
Moreover, the lower frame of Figs. 3 shows that, due to the non-negligible threshold energy the 
conventional method is (much) more unsensitive for lighter WIMPs (in contrast, the uncertainty 
interval given by Eqs. (17) and (18) becomes only a bit wider) and the curves thus go sharply 
upwards as the WIMP mass decreases. The incompatibility between two analyses becomes larger 
and one could therefore even give more strict constraints on the WIMP mass and the SI cross 
section (m^ > 45 GeV and > 7.5 x 10~^ pb in our simulation). 

In Figs. 4 we examine the same comparison of two analyses for a rather light input WIMP 
mass of 25 GeV. The lower frame shows that, with the non-negligible threshold energy one 
could even give the upper bounds of the WIMP mass and its cross section on protons (nucleons) 
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Figure 3: The SI WIMP-proton cross section cr^p estimated by Eqs. (17) and (12) and the lower 
and upper bounds of its Icr statistical uncertainty as functions of the WIMP mass (dashed blue 
curves) for a ^^Ge target. The input WIMP mass is 100 GeV. The threshold energies have been 
set as (upper) and 5 (lower) keV, respectively. The four extra curves have been drawn conven- 
tionally by using the shifted Maxwellian velocity distribution with four Sun's orbital velocities: 
fo = 180 km/s (dash-double-dotted orange), vq = 200 km/s (dash-dotted cyan), vo = 220 km/s 
(dotted magenta), vo = 240 km/s (double-dotted black), and the form factor given in Eq. (20). 
The other parameters are as in Fig. 2. See the text for further details. 
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Figure 4: As in Figs. 3, except that the input WIMP mass is only 25 GeV here. 



(20 GeV < < 50 GeV and 7 x lO"*^ pb < af^ < 1.1 x 10"^ pb in our simulation) from 
a single experiment. However, remind that with a fixed maximal cut-off energy and a number 
of total events, the higher the threshold energy, the larger the required exposure. Moreover, 
as we can see in the lower frame of Figs. 4, the prefactor, or equivalently, Iq, in Eq. (17) is 
underestimated due to a low (in contrast to the case shown in Figs. 1 and 2) kinematic maximum 
of the recoil energy. For a WIMP mass of 25 GeV, this kinematic maximum for a Ge target 
is 52.6 keV. Remind also that the recoil energy spectrum is approximately exponential, thus 
between Q = and Q = Qmax.kin = 52.6 keV, only ~ 53% of the total events are with energies 
Q > Qmin = 5 keV. Due to this underestimate of Iq, we will see later that the non-negligible 
threshold energy could cause serious problem once the WIMPs are (pretty) light. 

3 Estimating the SI WIMP— nucleon coupling 

In this section I consider further the case that two (or more) experiments with different target 
nuclei observe positive WIMP signals. 

3.1 Combining different experiments 

In Figs. 5 I show the SI WIMP-proton cross sections o'^p(m^) estimated by Eqs. (17) and (12) as 
functions of the WIMP mass on the a^^ — m^ plane for four different target nuclei: ''^Ge (dashed 
blue), ^^Si (dotted magenta), ^°Ar (dash-double-dotted orange), and ^^^Xe (long-dash-dotted 
cyan). Not surprisingly, all four curves pass through (approximately) the same values of 
and cr^p. It is in fact one of the basic ideas of the model-independent determination of the 
WIMP mass [6, 7] mentioned in the introduction^. However, one can also find here that the 
(approximately) common values of and cr^p are somehow underestimated, especially for the 
heavier input WIMP mass (see the lower frame). In Ref. [7], we discussed this phenomenon and 
introduced therefore an algorithmic procedure to correct this systematic deviation by matching 
the maximal cut-off energies of different targets. In Figs. 5 the vertical double-dashed black fines 
show the reconstructed WIMP masses and the lower and upper bounds of their la statistical 
uncertainties estimated by this algorithmic procedure. 

Once the WIMP mass on the right-hand side of Eq. (17) can be determined by means of 
the model-independent method with two different target nuclei, one can estimate the SI WIMP- 
proton coupling (cross section) straightforwardly. Here r{Qram) and Iq in the prefactor can be 
estimated from either one of the data sets used for determining or a third (independent) 
experiment. In Figs. 6, I show the reconstructed spin-independent WIMP-proton coupling 
I/pIj^pj, as a function of the input WIMP mass m^^i^. The expected number of total events in 
each data set has been set as 50 events on average under the experimental maximal cut-off 
energy Qmax set as 100 GeV for all targets; the experimental threshold energies are assumed to 
be negligible. Four nuclei: ^^Ge, ^^Si, ^°Ar, and ^^^Xe have been chosen for estimating r((5min) 
and Jo in Eq. (17). Following our work on the determination of the WIMP mass [7], ^^Si and 
'''^Ge have been chosen as two target nuclei for estimating m^, in Eq. (17). 

As a comparison to the use of the reconstructed WIMP mass (solid red), we consider here 
also the case that the WIMP mass m^. in Eq. (17) can be determined from some other (collider) 
experiments (dashed blue) with a higher precision. The input (true) WIMP mass has been used 
with an overall uncertainty of 5% for this case. Note that, firstly, in order to avoid complicated 
calculations of the correlations between the uncertainty on estimated by the algorithmic 



'A brief review of the determination of the WIMP mass is given in the appendix. 
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Figure 5: The SI WIMP-proton cross sections cT^p(^x) estimated by Eqs. (17) and (12) as 
functions of the WIMP mass for four different target nuclei: ™Ge (dashed blue), ^^Si (dotted 
magenta), "^^Ar (dash-double-dotted orange), and ^^^Xe (long-dash-dotted cyan). The input 
WIMP mass has been set as 100 (upper) and 300 (lower) GeV. The threshold energies for all 
targets are assumed to be negligible. The vertical double-dashed black lines show the recon- 
structed WIMP masses and the lower and upper bounds of their Icr statistical uncertainties 
estimated by the algorithmic procedure introduced in Ref. [7]. The other parameters are as in 
Fig. 2. 




-"Si (+ ^^si + ^'^ge), Q^3, < 1 00 keV, 1 (3) x 50 events, c^p^' = 1 o" pb 



AMIDAS http://pisrvO.pitphysik.Lni-tLebingen.de/darkmatler/amidas/ 




50 70 100 



500 700 1000 



50 70 100 



500 700 1000 



'"Ar (+ ^°si + '"ge), Q^a< < 1 00 keV, 1 (3) x 50 events, o^, 

AWIDAS http://pisrv0.pit.physik.Lni-lLebingen.de.'darkmatter/amidas/ 



= 10 ° pb 




'Xe (+ '°si + '"ge), Q^^^ < 100 keV, 1 (3) x 50 events, <j„ 

AWIDAS http://pisrv0.pit.physik.Lni-lLebingen.de./darkmatter/amidas/ 



too 200 

i„ [GeV] 




too 200 

i„ [GeV] 



Figure 6: The reconstructed SI WIMP-proton couplings l/pl^ec and the lower and upper bounds 
of their la statistical uncertainties as functions of the input WIMP mass ?^X:in- The long-dashed 
green curve indicates the theoretical value of the SI coupling. The solid red and dashed blue 
curves indicate the reconstructed SI couplings estimated with the reconstructed and the input 
(with an overall uncertainty of 5%) WIMP masses. "^^Ge, ^^Si, ^'^Ar, and ^^^Xe four nuclei 
have been chosen for estimating r((5min) and Jq in Eq. (17). Following our work in Ref. [7], 
^®Si and ^^Ge (labeled in the plots with small letters for their chemical symbols to indicate 
the independence of these data sets of the data set of the first nucleus) have been chosen as 
two target nuclei for reconstructing the WIMP mass m^. Parameters are as in Figs. 5. Note 
that for m^ in > 500 GeV the upper bounds of the statistical uncertainty are systematically 
underestimated. See the text for further details. 



procedure and those on r(Qinin) and Iq, we have assumed here that the two data sets with the Ge 
(Si) nucleus are independent of each other^. Secondly, an upper cut-off limit on the reconstructed 
WIMP mass has been set as 3000 GcV in our simulation. But, due to the very few number of 
events, the upper bounds of the la statistical uncertainty on the reconstructed mass for heavier 
input masses excess this limit. Hence, in the heavy mass range (^x,in > 500 GeV) in Figs. 6 
(and also in Figs. 7) the upper bounds of the la statistical uncertainty on the reconstructed SI 
couplings with all four targets arc systematically underestimated. 

It can however be found in Figs. 6 that, firstly, the reconstructed coupling |/p|rec estimated 
with the input (true) WIMP mass (dashed blue curves) for all four targets are underestimated for 
WIMP masses > 100 GeV; for the heavier target nuclei, Ge and Xe, this deviation is larger 
than for the lighter nuclei. Si and Ar. This is caused by the underestimate of Iq in Eq. (17), 
which we found in Figs. 1 and 2 and discussed there. Secondly, due to an underestimate of the 
reconstructed WIMP mass^'', the reconstructed couplings |/p|rcc with the reconstructed WIMP 
mass (solid red curves) for all four targets are more strongly underestimated than those with 
the true WIMP mass for WIMP masses > 100 GeV. Moreover, for lighter WIMP masses 
("T.x,in ^ 100 GeV), the reconstructed |/p|rec all targets are also a bit underestimated. This 
could possibly be caused by the statistical fluctuation due to the pretty few events. 

The systematic deviation of the reconstructed couplings with both reconstructed and real 
WIMP mass caused by the underestimate of Iq reflects the fact that the heavier the target 
nucleus, the more the contribution from WIMPs with higher velocities to the recoil spectrum. 
This observation implies that lighter nuclei could be better for estimating Iq, and in turn for 
reconstructing the SI WIMP-nucleon coupling. This can be seen more clearly from the difference 
of the reconstructed |/pp with Si, Ar and Ge, Xe for the case with the input (true) WIMP mass 
(blue dashed curves). However, Figs. 6 show also that the statistical uncertainties on |/pp 
estimated with the fighter nuclei are a bit larger than those with the heavier nuclei. And for 
heavier WIMP masses the deviation of Iq could in principle be alleviated by extending the 
experimental maximal cut-off energy Qmax to higher energy ranges, as discussed in the previous 
section. Moreover, remind that we simulated here with the same expected event number for all 
four target nuclei. In practice, we could measure (much) less WIMP events in experiments with 
fighter target nuclei {dR/dQ oc A^). This indicates also a larger statistical uncertainty. 

Nevertheless, our simulations shown in Figs. 6 demonstrate that, firstly, in spite of the 
systematic deviation for heavier WIMP masses due to the underestimate of /q, the true value of 
l/pp always lies within the la statistical uncertainty intervals. Secondly, for a WIMP mass of 
100 GeV, one could in principle estimate the squared SI WIMP-proton coupling with a statistical 
uncertainty of ~ 40% for the Si and Ar targets or of only ~ 30% for the Ge and Xe targets 
with only 50 events from one experiments^. This is much smaller than the uncertainty on the 
estimate of the local Dark Matter density (of a factor of 2 or even larger) . 



^The formulae needed for calculating the correlations between the uncertainties on the prefactor and on the 
WIMP mass estimated by two basic expressions given in Eqs. (31) and (A21) {not by the algorithmic procedure) 
are given in the appendix. 

^°The WIMP mass has been recostructed by a different program than that used in Ref. [7]. 

^^Note that these uncertainties have been estimated by 



la upper/lower bound of |/pp — |/p 
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rec 



(22) 



Pirec 



and, as shown in Figs. 6 as well as in Figs. 7, are asymmetric since the upper/lower uncertainties on the 
reconstructed WIMP mass are asymmetric [7]. 
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Figure 7: The reconstructed SI WIMP-proton couplings |/p|rec and the reconstructed WIMP 
mass m^^rec estimated by the method described in Ref. [7] with the Si and Ge targets on the 
cross section (coupling) versus WIMP mass plane. The filled green squares indicate the input 
WIMP masses and the theoretical values of the SI coupling. The red crosses (filled blue circles) 
indicate the reconstructed (input) WIMP masses and the reconstructed SI couplings estimated 
with these WIMP masses. The horizontal (vertical) solid red and dashed blue lines show the Icr 
statistical uncertainties on rn^^j-cc (l/plrec)- Parameters are as in Fig. 6. 



3.2 Constraints on the cross section— mass plane 



Since by using two (or three) data sets the squared SI WIMP-nucleon couphng |/pp can be 
estimated from experimental data directly without knowing the true value of the WIMP mass 
m^, and the same data sets can also be used to reconstruct m^, one can practically combine 
the reconstructed |/pp with the reconstructed together on the cross section-mass plane. In 
Figs. 7 I show the reconstructed SI coupling |/p|rg(, and the reconstructed WIMP mass m^^rec 
estimated by the algorithmic procedure described in Ref. [7] with the Si and Ge targets on 
the cross section (couphng) versus WIMP mass plane. The horizontal (vertical) sohd red and 
long-dashed bhie hues show the la statistical uncertainties on m^^^ec (l/plrec)- can be seen 
that the la statistical uncertainty areas of the reconstructed WIMP mass and its coupling can 
always cover their true values up to an input mass of ~ 1 TeV, although both of them are 
underestimated . 

The emphasis here is that, while by the conventional analyses for determining the WIMP 
mass and its SI coupling on nucleons (see e.g., [21, 22, 23]) one needs a model of the velocity 
distribution of halo WIMPs, one can estimate and |/pp separately by the method presented 
here with neither prior knowledge of each other nor an assumption about the WIMP velocity 
distribution. Certainly, how well one can estimate these two quantities depends not only on the 
event number but also on the target nucleus, as discussed in Ref. [7] and shown in Figs. 6 and 7. 

In Figs. 3 and 4 we saw that the non-negligible threshold energy could allow us to give 
more strict constraints on the WIMP mass and its SI coupling on nucleons. In Figs. 8 we 
therefore take into account a minimal cut-off energy Qmm = 5 keV for the first Ge target 
used for estimating r((5min) and Iq. It can be seen obviously that, for lighter WIMP masses 
("^x ~ 50 GeV) the reconstructed coupling |/p|rec is (strongly) underestimated for both cases 
with the reconstructed and the input (true) WIMP masses. As discussed at the end of the 
previous section, this is caused by a (very) low kinematic maximum of the recoil energy and, 
consequently, the underestimate of Iq. For a WIMP mass of 10 GeV, this kinematic maximum is 
just 11.8 keV and between Q = and Q = Qmax.km = 11-8 keV, only ~ 6.4% of the total events 
are with energies Q > Qmin = 5 keV! In contrast, for heavier WIMP masses (m^ > 50 GeV), 
the non-negligible threshold energy causes only slightly larger statistical uncertainties on the 
reconstructed SI couplings. 

So far we have assumed that each experiment "only" has an exposiire corresponding to 50 
total events. In Figs. 9 we raise this number by a factor of 10. Not surprisingly, all uncertainties 
on both the reconstructed WIMP mass and the reconstructed SI couplings shrink by a factor 
> 3 compared to the results shown in Figs. 7. Moreover, the small underestimate for lighter 
WIMP masses (m^ < 50 GeV) found in our simulations with only 50 events (see Figs. 6) 
disappears now. Note that, for heavier WIMP masses (m^ > 500 GeV), the upper bounds 
of the la statistical uncertainty on the reconstructed WIMP masses is now down to below our 
cut-off limit. 



4 Estimating the SD WIMP— nucleon couplings 



For the sake of completeness, I consider briefly in this section the case that the spin-dependent 
WIMP-nucleus interaction dominates over the spin-independent one. Then the WIMP-nucleus 
cross section ao in Eq. (2) can be expressed as [1, 2]: 
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Figure 8: As in Figs. 6 and 7, except that the experimental minimal cut-off energy for the first 
Ge target has been set as 5 keV. 




Figure 9: As in Figs. 6 and 7, except that the expected number of total events in all three 
experiments has been set as 500. 



Here Gp is the Fermi constant, J is the total spin of the target nucleus, (-S'(p,n)) are the ex- 
pectation values of the proton and neutron group spins, and a(p,n) are the effective SD WIMP 
couplings on protons and on neutrons. For the SD WIMP-nucleus cross section, it is usually 
assumed that only unpaired nucleons contribute significantly to the total cross section, as the 
spins of the nucleons in a nucleus are systematically anti-aligned^^. Under this assumption, the 
SD WIMP-nucleus cross section given above can be reduced to 
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Since for a proton or a neutron J — ^ and (^p) or (^n) 
protons or on neutrons can be given as 



|, the SD WIMP cross section on 
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By comparing Eq. (24) with the second expression in Eq. (11), the squared SD WIMP couplings 
on protons and on neutrons can be obtained from Eq. (17) straightforwardly as 
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Note that, for estimating Iq here by using Eq. (15), the elastic nuclear form factor F^(Q) must 
be chosen for the SD interaction. 

As the use of Eq. (17) for constraining the SI coupling |/pp discussed in Sec. 2, by assuming 
that the SD WIMP-nucleus interaction contributes to the total cross section ao dominantly, and 
combining with the conventional cross section-WIMP mass analysis, one could in principle use 
Eq. (26) to give the lower bounds of the WIMP mass and its SD couplings on protons and on 
neutrons from a single experiment with a target nucleus having spin sensitivity (almost) only on 
protons or on neutrons. By comparing these results with those obtained from the SI case, one 
could examine backwards the assumption for a dominant SD WIMP interaction. Meanwhile, as 
discussed in Sec. 3, once the WIMP mass ruy^ can be determined, one could then estimate the 
SD WIMP-nucleon cross sections by Eqs. (26) and (25) straightforwardly. 

Furthermore, by combining two target nuclei, one (X) of them has (almost) only spin sensi- 
tivity on protons and the other one (Y) on neutrons, one can easily find an expression for the 
ratio between two SD WIMP-nucleon couphngs from Eq. (26) as 
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^^However, more detailed nuclear spin structure calculations show that the even group of nucleons has some- 
times also non-negligible spin [1]. 
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'R-a,Y, T^j,n,Y, and TZn,Y can be defined analogously. Here m(x,Y) F{x,y){Q) are the masses 
and tfie form factors of the nucleus X and y, respectively, r(x,y)(Qmin,(x,y)) refer to the counting 
rates for the target X and Y at the respective lowest recoil energies included in the analysis, 
and £{x,Y) are the experimental exposures with the target X and Y . For the cancellation of the 
factors involving rriy. in the first line of Eq. (27), I used the general estimator for the WIMP 
mass given in Refs. [6, 7]: 
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Note that Eq. (27) derived here is in fact a special case of the general expression (16) given 
in Refs. [18, 17]. Detailed discussions about model-independent determinations of the ratios 
between different WIMP-nucleon couphngs/cross sections can be found in Refs. [18, 17, 19]. 



5 Summary and conclusions 

In this paper I presented the method for estimating the spin-independent WIMP-nucleon cou- 
pling from elastic WIMP-nucleus scattering experiments. This method is independent of the 
velocity distribution of halo WIMPs as well as (practically) of the as yet unknown WIMP mass. 
Assuming that an exponential-like shape of the recoil spectrum is confirmed from experimental 
data, the required information are only the measured recoil energies and the number of events 
in the first energy bin from at least two experiments with different target nuclei as well as the 
unique assumption for the local WIMP density. 

In Sec. 2 I rederived the expression for estimating the (squared) SI WIMP-nucleon coupling 
l/pP as a function of the (unknown) WIMP mass [17]. Then I demonstrated that, by comparing 
the constrained area estimated by this method to that given by the conventional analysis with 
an assumed halo model, one could in principle - for the first step with only one experiment 
observing positive signals - give the lower bounds of the WIMP mass and its SI cross section on 
nucleons from a single experiment. 

For the next step, I discussed in Sec. 3 that, by using measured recoil energies from two (or 
three) experiments with different target nuclei, we could not only determine the WIMP mass as 
discussed in Refs. [6, 7], but also estimate the SI WIMP-nucleon coupling, with neither prior 
knowledge of each other nor an assumption for the velocity distribution of halo WIMPs. 

However, due to the degeneracy between the local WIMP density and the WIMP-nucleus 
cross section, it is impossible to determine both of them independently. As the simplest way 
one has thus to make an assumption for the local WIMP density. Nevertheless, since the SI 
WIMP-nucleon coupling is inversely proportional to the local WIMP density, whose common 
value would possibly be underestimated, one can then at least give an upper bound on this 
couphng. Moreover, our simulations show that, in spite of the very few (0(50)) total events 



from one experiment, for a WIMP mass of 100 GeV, the SI WIMP-nucleon coupling can be 
estimated with a statistical uncertainty of only ~ 15%; it leads to an uncertainty on the SI 
WIMP-nucleon cross section of only ~ 30%, which is (much) smaller than the uncertainty on 
the estimate of the local Dark Matter density (of a factor of 2 or even larger). 

Our simulations show also that, due to (mainly) the experimental maximal cut-off energy, 
the SI WIMP couphng could be underestimated for heavier WIMP masses, especially with 
heavy target nuclei, e.g., Ge or Xc. However, since the kinematic maximum of recoil energies 
for heavier WIMP masses and/or with heavy target nuclei are (much) higher than for lighter 
WIMP masses with light nuclei, one could practically alleviate this systematic deviation by 
extending the detector sensitivity to higher energy ranges. Moreover, due to the fairly large 
statistical uncertainty, the true value of the SI WIMP-nucleon couphng lies always within the 
la statistical uncertainty interval. 

In Sec. 4 1 turned to consider the case that the spin-dependent WlMP-nucleus interaction 
dominates over the SI one. By assuming (naively) that only unpaired nucleons contribute sig- 
nificantly to the total WIMP-nucleus cross section, I gave also the expression for estimating the 
(squared) SD WIMP-nucleon couplings |a(p,n)P as functions of the (unknown) WIMP mass. As 
for the SI case, by comparing the constraints estimated by this method to those given by the 
conventional analysis, we could in principle also give the lower bounds of the WIMP mass and 
its SD cross sections on nucleons from a single experiment. 

Our simulations presented here are based on several simplified assumptions. Firstly, the 
sample to be analyzed contains only signal events, i.e., is free of background^^' Secondly, 
all experimental systematic uncertainties as well as the uncertainty on the measurement of the 
recoil energy have been ignored. The energy resolution of most currently running and projected 
detectors is so good that its uncertainty can be neglected compared to the statistical uncertainty 
with (very) few events in the foreseeable future. 

A non-negligible threshold energy makes the conventional model-dependent analysis less 
sensitive on light WIMPs {rriy- < 20 GeV), it could however give us more strict constraints on the 
lower bounds of the WIMP mass and its couplings on nucleons. In contrast, our simulation shows 
that, by using our model-independent method, the non-negligible threshold energy could cause 
not only a larger statistical uncertainty on the reconstructed couplings, but also a significant 
underestimate if WIMPs are (very) hght. 

In summary, I demonstrated in this paper the use of our new method for estimating the 
spin-independent WIMP-nucleon coupling with neither a prior knowledge of the WIMP mass 
nor an assumption for the velocity distribution of halo WIMPs. By combining with information 
on the ratios between different WIMP-nucleon couplings/cross sections, which could also be 
determined model-independently [18, 17, 19], one could in principle also estimate the absolute 
values of the spin-dependent cross sections. This information combined with the reconstructed 
WIMP mass will allow us not only to constrain the parameter space in different extensions of 
the Standard Model of particle physics [30, 31, 16], but also to identify WIMPs among new 
particles produced at colliders [8]. Furthermore, knowledge of the WIMP mass and its couplings 
could not only offer a new approach for estimating the local WIMP density, but also permit the 
prediction of the WIMP annihilation cross section and the event rate in indirect Dark Matter 
detection experiments [1, 2]. 

^^For background discrimination techniques and status in currently running and projected direct detection 
experiments see e.g., [24, 25, 26, 27]. 

^^For detailed simulations and discussions about effects of residue background events on the reconstructions of 
the WIMP mass and its SI coupling on nucleons see [28, 29]. 
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A Formulae needed in Sees. 2 and 3 

Here I list all formulae needed for our model-independent data analyses described in this article. 
Detailed derivations and discussions can be found in Refs. [5, 7]. 

A.l Estimating r(Qmin) and /n(Qmin, Qmax) 

First, consider experimental data described by 

Qn-k<Qn,i<Qn+k^ ^ = 1, 2, • • • , iV„, n = 1, 2, • ■ ■ , S. (Al) 

Here the total energy range between Qmin and Qmax has been divided into B bins with central 
points Qn and widths 6„. In each bin, Nn events will be recorded. Since the recoil spectrum 
dR/dQ is expected to be approximately exponential, the following ansatz for the measured recoil 
spectrum ( before normalized by the experimental exposure £) in the nth bin has been introduced 
[5]: 

/dR\ ^(^\ = r„ e^"^'^-^^'"^ (A2) 

Here r„ is the standard estimator for {dR/dQ)expt at Q = Qn- 

rn = ^, (A3) 

kn is the logarithmic slope of the recoil spectrum in the nth Q— bin, which can be computed 
numerically from the average value of the measured recoil energies in this bin: 

= (I) coth (^A - i- , (A4) 
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The error on the logarithmic slope kn can be estimated from Eq. (A4) directly as 
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with 



Qs,n in the ansatz (A2) is the shifted point at which the leading systematic error due to the 
ansatz is minimal [5], 



Qs,n Qn I, 



(A8) 



Note that Qs,n differs from the central point of the nth bin, Q^. Prom the ansatz (A2), the 
counting rate at Q — Qmm can be calculated by 



and its statistical error can be expressed as 
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Finally, since all are determined from the same data, they are correlated with 

iVtot Q{n+m~2)/2 

cov(4,/J = E^^, 
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(A12) 



where the sum runs over all events with recoil energy between Qmin and Qmax- And the correlation 
between the errors on r{Q-[am), which is calculated entirely from the events in the first bin, and 
on /„ is given by 
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note that the sums here only count in the first bin, which ends at Q = Qmin + ^i- 

On the other hand, with a functional form of the recoil spectrum (e.g., fitted to experimental 
data), {dR/ dQ)e:x.pt, one can use the following integral forms to replace the summations given 
above. Firstly, the average Q— value in the nth bin defined in Eq. (A5) can be calculated by 
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For /n(Qmm, Qmax) givcn in Eq. (15), we have 

/•Qmax g("-i)/2 /dR\ 

-^n(Qmin) Qmax) — / ,-,9 dQ , 

and similarly for the covariance matrix for 7„ in Eq. (A12), 
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Remind that (di?/d(5)expt is the measured recoil spectrum before normahzed by the exposure. 
Finally, /^(Qmin, <5min + ^'i) needed in Eq. (A13) can be calculated by 
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(A17) 



Note that, firstly, r((5min) and /„((5min, Qmin + W) should be estimated by Eqs. (A9) and (A17) 
with ri, ki and Qs,i estimated by Eqs. (A3), (A4), and (A8) in order to use the other formulae for 
estimating the (correlations between the) statistical errors without any modification. Secondly, 
'"(Qmin) and /n(Qmin, Qmax) estimated from a scattering spectrum fitted to experimental data 
are usually not model-independent any more. Moreover, for estimating the SD WlMP-nuclcon 
couplings by Eq. (26), the elastic nuclear form factor F^{Q) in Eqs. (15), (A12), (A15), (A16), 
and (A17) should be understood to be chosen for the SD interaction. 



A. 2 Determining the WIMP mass m 



X 



By requiring that the values of a given moment of fi{v) estimated by Eq. (13) from two exper- 
iments with different target nuclei, X and Y, agree, appearing in the prefactor a" on the 
right-hand side of Eq. (13) can be solved analytically as [6, 7]: 
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with TZn,{x,Y) given by Eq. (30). Note that the general expression (31) can be used either 
for spin-independent or for spin-dependent scattering, one only needs to choose different form 
factors under different assumptions; the form factors needed for estimating In,{x,Y) by Eq. (15) 
or (A15) are thus also different. 

By using the standard Gaussian error propagation, a lengthy expression for the statistical 
uncertainty on m^\^^„^ can be obtained as 
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Here a short-hand notation for the six quantities on which the estimate of rriy^ depends has been 
introduced: 
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and similarly for the Q^y. Estimators for cov(ci,Cj) have been given in Eqs. (A12) and (A13). 
Explicit expressions for the derivatives of TZn,x with respect to Ci^x are: 



dTZn,x _ n + 1 
dln,x n 

dnn,x ^ _i 

dIo,x n 



min,X J 



2Q!:in'i^Vx(gmin,x) + {u + l)/n,xFl (Qmin.x) 



7^, 



n,X 



1 /2 2 
2<5min,x'^x(<5min,x) + -^0,X-^l(<5min,x) 



7^n 



X 



(A20a) 



(A20b) 



and 



■n,X 



drx{Qmin,x) 



2 

n 



2Qt^,Trx{Qmu.,x) + (n + l)In,xFUQmin,x), 



X 



,mm,X I 



2QmiiiX^^(Qmin,x) + Io,X Fx{Qmm,x) 



n,X 



(A20c) 



explicit expressions for the derivatives of 7^„,y with respect to Cj^y can be given analogously. Note 
that, firstly, factors 7ln,{x,Y) appear in all these expressions, which can practically be cancelled 
by the prefactors in the bracket in Eq. (A18). Secondly, all the Io,(x,y) ^n,(x,Y) should be 
understood to be computed according to Eq. (15) or (A15) with integration limits Qmin and 
Qmax specific for that target. 

On the other hand, since |/pp in Eq. (17) is identical for different targets, it leads to a second 
expression for determining [7]: 

_ {mx/myf^'^ my - mxiTl^,x/T^a,Y) 



Tla,x/'^a,Y - {rnx/mYf''^ 



(A21) 



Here oc ^(x,y) has been assumed, and TZa,{x,Y) have been given in Eq. (28). Similar to 

the analogy between Eqs. (31) and (A21), the statistical uncertainty on m^|^ can be expressed 
as 



{mx/myf/' 


ITLX — fny 


(7^,,x/7^^,y) 




Tia,x/Tla,Y - (mx/mYf^^ 


2 



X 



E 



x 



'<-a,X i,j=2 



dc 



LX 



cov(ci,x, Cj^x) + {X — >Y) 



1/2 



(A22) 



where I have used again the short-hand notation in Eq. (A19); note that Ci^{x,y) — In,{x,Y) do 
not appear here. Expressions for the derivatives of TZ^j^x can be computed from Eq. (28) as 

dTZ^^x _ FxiQm.m,x) 



9Io,X [2g^f„_^rx((5mm,x) + Io,xF^{Qram,x) 



n. 



a,X , 



on. 



cr,X 



drxiQ 



min,X J 



2Qmi 



1/2 



m'm.X 



1 /2 

2(3min,X^x((3min,x) + Io,xFx{Qmin,x) 



a,X 



(A23a) 
(A23b) 



and similarly for the derivatives of TZa,Y- Remind that factors TZa,{x,Y) appearing here can also 
be cancelled by the prefactors in the bracket in Eq. (A22). 

In order to yield the best-fit WIMP mass as well as to minimize its statistical uncertainty 
by combining the estimators for different n in Eq. (31) with each other and with the estimator 
in Eq. (A21), a function has been introduced as [7] 



X'(rnx) = E ikx - Ay) C-^ (hx - Uy) , 



(A24) 



where 



i,X 



a 



X 



2gSx'^x(Qmin)/Fl(gmin,x) + (i + l)/.,X 



2QX,X^x(gmin)/Fl(Qmin,x) + Io,X 

oixT^ix y 



300 km/s 



300 km/s^ 



(A25a) 



for i = -1, 1, 2, . . . , rimax, and 



X 



42 



2Qmin, x'^x(Qmin)/-f'x(Qmin,x) + /o,X 



42 



T^a,x \m^ + mx. 



(A25b) 



the other nmax + 2 functions /i,y can be defined analogously. Here nmax determines the highest 
moment of fi{v) that is included in the fit. The are normalized such that they are dimen- 
sionless and very roughly of order unity in order to alleviate numerical problems associated with 
the inversion of their covariance matrix. Note that the first Timax + 1 fit functions depend on 
only through the overall factor a and in Eqs. (A25a) and (A25b) is now a fit parameter, 
which may differ from the true value of the WIMP mass. Finally, C in Eq. (A24) is the total 
covariance matrix. Since the X and Y quantities are statistically completely independent, C can 
be written as a sum of two terms: 



dj = COV {fi^x, fj,x) + GOV (/i,y, fj^y) 



(A26) 



The entries of the C matrix given here involving basically only the moments of the WIMP 
velocity distribution can be read off Eq. (82) of Ref. [5], with an slight modification due to the 
normalization factor in Eq. (A25a)^^: 

COV (/,, /,) = [/, cov(7o, lo) + + W + l)cov(7,, /,) 

- a^{j + l)fi cov(/o, Ij) - a\i + l)fj cov(/o, h) 

+ DiD^a\r{Q^^)) - {DJj + D^h) cov(r(Q^i„), /q) 

+ &{j + 1) A cov(r(g^in), Ij) + a\i + l)Dj cov{r{Q 



mm J y -^i] 



Here I used 



1 



2QLfnKQmin)/i^^(Qmin)+/0 

a 



a = 



and 



D, 



300 km/s ' 

^f^ 



1 



dr(Qmm) 

. . . , Tim; 

2 



F (Qmin) 

and 



Q! Qrai 



(i+l)/2 
min 



for i = -1, 1, 2, 

-^nmax+l ~ "E^tTT^ \ (~Qmin/) 



'T'ma.x~hl 



)■ 



(A27) 

(19) 
(A28) 

(A29a) 
(A29b) 



Finally, since the basic requirement of the expressions for determining given in Eqs. (31) 
and (A21) is that, from two experiments with different target nuclei, the values of a given 



^^Since the last fi defined in Eq. (A25b) can be computed from the same basic quantities, i.e., the counting 
and the integrals Iq, it can directly be included in the covariance matrix. 



moment of the WIMP velocity distribution estimated by Eq. (13) should agree, the upper cuts 
on fi{v) in two data sets should be (approximately) equal^^. Since Vcut — OiVQmax, it requires 
that [7] 



Q 



max,y 



Kay J 



2 

) Qmax,X ■ 



(A30) 



Note that a defined in Eq. (5) is a function of the true WIMP mass. Thus this relation for 
matching optimal cut-off energies can be used only if is already known. One possibility to 
overcome this problem is to fix the cut-off energy of the experiment with the heavier target, 
minimize the x^{it^x) function defined in Eq. (A24), and then estimate the cut-off energy for 
the fighter nucleus by Eq. (A30) algorithmically [7] . 



A. 3 Covariance of and l/ATm 

First, the statistical error on l/A/'m can be given from Eq. (19) directly as 



a 



-^^(Qmin) 



-^^(Qmin) 



cov(r(g^in),/o)- (A31) 



For the case that one has only two data sets with different target nuclei, X and F, one of these 
two data sets will then be needed for reconstructing the WIMP mass and also for estimating 
l/A/'m in Eq- (17). The uncertainties on and l/A/'m are thus correlated. Assuming that the 
WIMP mass is reconstructed by Eq. (31), and target X{Y) is used for estimating l/A/'m, the 
covariance of ^xl(i)"> 1/-Mn,(x,y) can be obtained by modifying Eq. (A18) slightly as 

cov(mx|^^„^,l/A/'m,x) 



('7^n,x/^n,y - \imx 
i=i \ 9ci,x 



3 



/my) 

COv(Ci,x, Io,x) + COv(q,x, rxiQmm,x)) 



2Qmin,X 



min,X I 



(A32a) 



and 



cov(m^|^^„^ , l/A/'m,y) 
^Jmx/my {mx - my) (7^n,x/^ 



(nn,x/T^n,Y - \lmx/my) 
dUnV 



3 
i=l 



dci 



i,Y 



COv(Ci,y, /o,y) + COv(Ci,y, Ty (Qmin.y)) 



min,y I 



(A32b) 



For the case that the WIMP mass is reconstructed by Eq. (A21), one can also modify Eq. (A22) 
to obtain that 

cov(m^|^,l/A/'m,x) 



Here the threshold energies have been assumed to be negUgible. 



{mx/niyY''' (nix - niy) (TZ^.x /'^as) 



\lla,xl'TlaX - {rnx/mYf^ 




dc. 



COv(q,x, h,x) + COv(Cj,x, rx(Qmin,x)) 



^Vmin,X 
Fx{Qmm,x) 



, (A33a) 



and 



{mx/mYf^'^{mx-mY){Tla,x/'^a,Y) / -1 



\lla,x/n^,Y - 



7^„ 




dc 



LY 



mx/mY)^^'^ 

COv(Ci,y, /o,y) + COv(Ci,y, ry((5inin,y)) 



20'/' 
^Vmin,y 



min,y J 



(A33b) 



Note that, firstly, in the above expressions we have to use {mx — ttiy) instead of |mx — ttiyI in 
Eqs. (A18) and (A22); for expressions with the Y target, there is an additional "— (minus)" 
sign. Secondly, the algorithmic process for matching the experimental maximal cut-off energies 
of two experiments used for the reconstruction of the WIMP mass can also be used with the 
basic expressions (31) and (A21). For this case and the hghter nucleus is used for estimating 
l/jVm, the energy range of the sum in Eq. (A12) or of the integral in Eq. (A16) as the estimator 
for the covariance of /„ should be modified to be between Qmin and the reduced maximal cut-off 
energy of the lighter nucleus. 
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